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By  computer  simulation  the  processes  of  porous  silicon  formation  and  of  silicon-on- 
porous  silicon  epitaxy  are  studied.  The  model  of  electrochemical  etching  is  applied  for 
p+Si  and  p'Si  substrates  and  takes  into  account  the  non-homogeneous  surface  charge 
distribution,  the  thermal  generation  of  holes  and  quantum  confinement  effects.  The 
epitaxy  is  studied  on  the  basis  of  Gilmer  model  extended  to  the  case  of  relief  surface. 
The  porous  structures  obtained  by  computer  simulation  are  presented  for  various  ano- 
dization conditions.  The  analysis  of  3D  images  shows  the  profile  of  porosity  over 
depth  and  multifractal  properties.  It  is  shown  that  the  mechanism  of  epitaxy  on  PS 
(111)  surface  is  provided  by  the  thin  pendant  layer  formation.  The  dependence  of  ki- 
netics of  epitaxy  upon  the  porosity,  molecular  flow  density  and  initial  surface  rough- 
ness are  established. 


1.  Introduction 

Porous  silicon  (PS),  prepared  by  anodization,  has  become  a very  promising 
material  for  device  application  due  to  its  visible  luminescence  [1],  extremely  large 
specific  surface  and  wide  technological  opportunities  [2,  3],  It  is  also  an  interesting 
subject  for  fundamental  investigation,  since  fractal  structure  of  PS  exhibits  a number 
of  optical,  transport  and  thermodynamic  properties  [4],  which  are  not  fully  understood. 

Along  with  experimental  methods,  the  computer  simulation  has  been  used  for 
the  study  of  PS  properties.  The  first  step  in  this  direction  was  the  modelling  of  PS 
formation.  The  diffusion  limited  aggregation  (DLA)  model  and  its  variations  [5,  6] 
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appeared  to  be  easily  adapted  for  computer  simulation.  However,  structures  obtained  in 
these  computational  experiments  could  hardly  be  used  in  further  investigation  mainly 
for  the  following  reasons.  First,  these  structures  were  modelled  on  the  2D  meshes. 
Second,  they  had  only  visual  resemblance  to  the  real  ones,  while  the  detailed  topologi- 
cal analysis  showed  that  much  of  structural  characteristics  (such  as  shape  and  size 
distribution  of  pores,  interface  profiles  etc.)  were  not  quite  adequate.  The  simulated 
porous  structures  obtained  by  Aleksandrov  and  Novikov  [7]  exhibited  the  variety  of 
morphologies  and  the  structural  [8]  and  thermodynamic  [9]  properties  of  PS  were 
studied.  The  qualitative  analysis  of  the  simulated  PS  layer  surface  was  used  for  model- 
ling of  silicon-on-porous  silicon  epitaxy.  The  present  paper  includes  the  series  of 
works  which  deal  with  the  modelling  of  PS  formation,  structural  analysis  of  various 
simulated  morphologies  and  homoepitaxy  on  PS. 

2.  The  model  and  method  of  calculation 

2.1.  THE  MODEL  OF  POROUS  SILICON  FORMATION 

Since  a hole  is  necessary  in  electrochemi- 
cal reaction  of  Si  dissolution  under  anodization, 
the  actual  question  is  at  what  points  a hole  pref- 
erably crosses  the  Si/HF  interface.  We  suppose 
that  the  morphology  of  the  porous  structure  is 
determined  by  self-consisting  correlation  be- 
tween potential  and  geometrical  relief  at  Si/HF 
interface.  At  the  first  stage  of  pore  formation  the 
initial  flat  and  smooth  silicon  anode  becomes  a 
relief  as  the  result  of  random  dissolution  (con- 
tour 1 in  Fig.  1).  This  stage  proceeds  until  tips  or 
cavities  at  the  Si/HF  interface  become  marked 
enough  to  trap  the  charge  carriers  (contour  2 in 
Fig.l).  F ions  are  trapped  in  the  cavities,  filled 
with  HF  solution,  while  holes  drift  to  the  hill- 
ocks (contour  3 in  Fig.l).  If  the  electrolyte  near 
the  silicon  anode  is  depleted  (high  current  den- 
sity approximation)  then  the  hillocks  are  the 
effective  centers  for  dissolution  and  electro- 
polishing occurs.  In  the  opposite  case  (low  current  density  and/or  light  doping  level 
and/or  high  HF  concentration)  electric  field  lines  are  focused  on  the  pore  tips,  where 
the  dissolution  preferably  takes  place.  In  an  intermediate  case,  when  electrolyte  is 
weakly  depleted,  only  some  randomly  formed  pores  induce  at  their  tips  an  electric  field 
large  enough  to  attract  the  holes  effectively.  The  current  density,  j,  at  the  tips  of  these 
pores  is  higher  than  the  average  over  the  surface  and  in  the  region  of  electrolyte  near 
the  pore  tips  the  strong  ion  depletion  is  formed.  This  leads  to  the  effect  of  “electro- 
polishing” in  micro  regions  near  pore  tips  rather  than  over  the  whole  surface.  The  size 


Fig.l.  Schematic  illustration  of  the 
pore  formation  relief  mechanism. 
The  contours  1-4  show  the 

consecutive  stages  of  Si/HF 

interface  evolution. 
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of  the  micro  regions  determines  an  average  diameter  of  pores,  which  is  constrained  by 
the  condition  of  the  depletion  in  the  electrolyte  at  the  pore  tip.  As  the  diameter  of  the 

pore  exceeds  some  critical  value  D,  the 
local  current  density  (under  the  fixed 
average  one)  becomes  too  low  to  provide 
the  strong  depletion  in  the  electrolyte. 
Then  the  region  of  “electropolishing” 
gets  narrow.  The  pore  diameter  value  is 
provided  by  the  self-consistent  processes. 

Computer  simulation  of  p+-Si  ano- 
dization is  performed  on  two-dimensional 
(160x160)  mesh.  Each  position  on  the 
mesh  is  defined  as  a cell.  One  mesh  step 
corresponds  to  the  real  size  of  surface 
charge  fluctuations.  The  computational 
algorithm  is  illustrated  in  Fig.2.  N holes 
undergo  random  transport  by  the  mesh, 
hopping  over  a distance  equal  to  the 
mean  free  path.  As  soon  as  a hole  is  lo- 
cated within  radius  RD  from  the  nearest 
pore  tip,  drift  in  the  direction  of  this  tip 
starts.  Rd  is  equal  to  the  Debye  radius  and 
expressed  in  mesh  cells.  After  the  hole  reaches  a point  on  the  interface,  the  analysis  of 
interface  curvature  is  carried  out  within  interval  L from  this  point.  The  dissolution 
occurs  in  the  most  marked  protuberance  over  the  interval.  The  average  diameter  of 
pores,  d,  was  estimated  roughly  as  follows  [7]: 


Fig.2.  The  illustration  of  the  computational 
algorithm  for  p+-PS  formation. 
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where  c (wt.%)  is  HF  concentration,  Ns  (cm'2)  is  the  surface  density  of  pores,  j 
(mA/cm2)  is  anodic  current  density,  T (K)  is  the  temperature,  EA= 343  meV,  C=3300 
mA  cm"2  wt.  %'3/2.  The  parameter  L (in  mesh  steps)  is  equal  to  d followed  by  Eq.  (2). 

In  lightly  doped  crystal  (p<1016  cm'3)  the  space  charge  density  and  the  local 
electrical  field  are  low.  Therefore,  the  effect  of  potential  relief  on  the  morphology  of 
porous  structure  is  expected  to  be  negligible  in  this  crystal.  For  this  reason  the  domi- 
nant mechanisms  of  the  pore  formation  are  diffusion,  thermal  generation  and  recom- 
bination of  holes  and  - in  nanocrystallites  - quantum  confinement  effect.  The  quantum 
confinement  is  known  to  increase  the  band  gap.  Consequently  a barrier  appears  be- 
tween the  bulk  region  and  nanoparticles,  which  prevents  hole  penetration  inside  nano- 
particles. The  size,  at  which  the  quantum  confinement  becomes  significant,  is  denoted 
Rq  and  equal  to  10  nm  and  less  [10], 

Computer  simulation  of  the  pore  formation  in  p'-Si  is  performed  on  two- 
dimensional  (160x160)  and  three-dimensional  (160x160x80)  meshes.  One  mesh  step 
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corresponds  to  the  average  interatomic  distance.  The  computational  algorithm  differs 
from  the  one  described  above.  Hole  may  disappear  from  the  current  point  and  appear 
in  another  position  of  crystal  volume  with  probability  G.  The  generation  rate  g (s'1)  is 
linked  to  G by  the  relationship: 

g=Gv</ 1,  (2) 

where  / ~10'6-10'7  m is  mean  free  path,  v0  ~105  m/s  - mean  velocity  of  hole.  After  a 
hole  reaches  the  interface,  the  analysis  of  surrounding  space  is  carried  out.  If  the  crys- 
tal volume  contains  the  sphere  of  radius  RQ,  then  the  dissolution  occurs. 

2.2.  THE  MODEL  OF  SILICON-ON-POROUS  SILICON  EPITAXY 

The  extended  diffusion  model  was  constructed  on  the  basis  of  the  initial  prem- 
ises of  the  Gilmer  model  and  includes  two  new  assumptions.  First,  vacancies  and  over- 
hangs are  permitted.  Second,  the  summation  of  the  second-nearest  neighbors,  which 
determine  the  activation  energy  of  diffusion,  extends  over  all  12  positions  in  the  second 
coordination  sphere.  Evaporation  of  atoms  is  neglected.  All  numerical  parameters  of 
the  extended  model  are  corrected  so  that  the  growth  process  on  the  smooth  Si  (111) 
surface  will  occur  identically  to  the  process  in  the  Gilmer  diffusion  model  [11]. 

The  computational  algorithm  is  performed  on  a A/t-xJVyx//z=160x  160x20  three- 
dimensional  grid  whose  nodes  correspond  to  regular  positions  of  atoms  in  the  crystal 
structure  of  silicon.  The  state  of  the  cells  that  are  filled  with  atoms  is  characterized  by 
the  numbers  Nj  and  N2  of  neighbors  in  the  first  and  second  coordination  spheres 
(l<//;<4,  1<A^12)  respectively.  Each  atom  executes  a diffusion  hop  with  probability  p 
(NltN2,T).  The  final  position  of  a diffusion  hop  is  chosen  randomly  among  the  unoccu- 
pied cells  in  the  first  and  second  coordination  spheres.  Transitions  which  result  in  an 
atom  not  having  any  neighbors  in  the  first  coordination  sphere  are  forbidden  (that  is 
why  Nj^O,  N2*Q).  The  probability  of  a diffusion  hop  (with  the  condition  that  the  end- 
point position  is  free)  equals 


P(N !,N2,T)=pcexp(-E/kT),  E=NrE1+NrE2,  (3) 

where  k is  Boltzmann's  constant,  T is  temperature,  and  p0=exp((E ,+3E2)lkl)  is 
a factor  chosen  so  that  the  probability  of  a diffusion  hop  on  the  free  (111)  surface 
equals  1.  To  decrease  the  amount  of  computation,  a threshold  probability  ps,  which 
limits  the  number  of  “candidates”  for  a diffusion  hop  (active  array)  to  atoms  for  which 
p(Ni,N2,T)>ps,  is  introduced.  Scanning  over  the  elements  in  an  active  array  is  per- 
formed during  iteration.  The  time  interval  corresponding  to  one  iteration  is  estimated 
as  v=%po/v,  where  v=2.5-1013  s'1  is  the  frequency  of  atomic  oscillations  and  %=6/16  is 
the  coefficient  that  takes  into  account  the  ratio  between  the  numbers  of  possible  final 
positions  for  Is(l  1 1)  in  the  Gilmer  and  the  extended  diffusion  models. 

The  parameter  E2  characterizes  the  kinetics  of  attachment/detachment  of  atoms 
on  the  faces  of  the  elementary  layers.  These  processes  were  investigated  by  simulation 
method,  and  the  value  E2= 0.2  eV  provided  good  agreement  with  experimental  data 
[12],  Since  in  our  model  an  adatom  on  a smooth  surface  has  one  neighbor  in  the  first 
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Rd=0 


Rd=10 


Rd=20 


Fig.  3.  A matrix  of  computer  simulation  plots  with  varying  effective  parameters  L and  Rd 
(iV=80). 

and  three  in  the  second  coordination  spheres,  the  activation  energy  of  surface  diffusion 
can  be  expressed  as  EA=E1+3E2.  At  the  same  time,  the  value  EA=  1.3±0.2  eV  was  ob- 
tained experimentally  in  a number  of  work.  We  employed  the  value  £^=1.34  eV,  which 
corresponds  to  E;=0.74  eV. 

The  probability  that  an  atom  from  a molecular  beam  reaches  the  surface  over 
time  x is 

Pif=2FxNxNy,  (4) 

where  F (bilayer/s)  is  the  molecular  flux  density. 

The  160x160x14  meshes  speckled  with  bottomless  cylindrical  pores  were  used 
as  the  initial  layer  for  simulation.  The  cylinders  were  randomly  distributed  over  the 
surface  and  oriented  normally  to  it.  The  number  and  radius  of  cylinders  were  the  pa- 
rameters of  simulation. 
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3.  Results  and  discussion. 

3.1.  RESULTS  OF  SIMULATION  OF  POROUS  SILICON  FORMATION. 

A matrix  of  plots,  obtained  by  computer  simulation  on  two-dimensional  mesh 
(160x160),  with  N=80  is  shown  at  Fig.3.  Black  colour  corresponds  to  pores.  The  pores 
grow  downwards  from  the  top  of  each  plot.  The  effective  parameters  RD  and  L (in  mesh 
steps)  are  shown  for  each  line  and  column  of  the  matrix.  The  pore  diameter  is  seen  to 

increase  with  the  increase  of  L, 
which  corresponds  to  current 
density  increase  and/or  HF 
concentration  decrease.  This 
result  agrees  with  experimental 
data  [13,  14],  Within  radius,  L, 
pores  are  continuos  in  plane- 
parallel  directions,  which  agrees 
with  high-resolution  and  cross- 
sectional  TEM  photographs  of 
the  mesoporous  Si  layers  [15], 
For  Z,=4  and  RD= 0 current  den- 
sity is  great  enough  to  achieve 
polishing.  RD  increase  leads  to 
straightening  of  the  pores  and 
increasing  distance  between 
them.  By  the  conditions  L= 0 
and  Rd=  1 0 and  20,  the  direction 
of  growth  is  practically  prede- 
termined, since  the  active  zone  for  each  pore  is  limited  by  its  lowest  point.  Similar 
conditions  were  provided  by  Lehmann  et.  al.  [16]  in  n-Si,  using  inducing  pits  and 
backside  illumination. 

The  three-dimensional  images  of  simulated  structures  (160x160x80)  were  ob- 
tained. The  calculations  show  that  the  porosity  P changes  over  the  depth,  h,  of  the 
porous  layer  in  accordance  with  SIMS  data  [17],  dP/dh  decreases  with  increasing  G. 
In  order  to  interpret  this  result,  it  is  necessary  to  first  explain  why  the  gradient  in  DLA 
model  rises  without  hole  generation  (G=0).  In  the  growth  process  limited  by  diffusion, 
holes  can  not  penetrate  from  the  bottom  to  the  top  of  the  growth  layer  since  they  are 
effectively  captured  by  the  branches  of  the  ramified  structure.  For  this  reason  the 
growth  preferably  proceeds  in  thin  front  layer.  Let  us  characterize  it  at  any  time  by  the 
current  active  surface  defined  as  the  number  of  possible  centres  of  growth  in  the  front 
layer.  The  initial  flat  Si/HF  interface  has  the  highest  active  surface.  As  the  cluster 
develops,  the  active  surface  decreases  (even  though  the  whole  specific  surface  may 
increase).  The  evolution  of  the  active  surface  determines  the  porosity  profile,  so  the 
change  in  active  surface  leads  to  the  porosity  gradient.  The  hole  generation  permits 


Structure  fragment  size  (mesh  cells) 


Fig.  4.  Fractal  dimension  of  structures  formed  by  com- 
puter simulation  on  three-dimensional  mesh  vs  size  of 
detail. 
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holes  to  appear  at  any  depth  of  PS  layer,  and  the  active  surface  is  no  longer  limited  by 
the  front  layer.  Then  the  pore  diameter  increases  until  the  quantum  confinement  in 
crystal  nanoparticles  becomes  essential.  So  the  simultaneous  action  of  hole  generation 
and  quantum  confinement  leads  to  smoothing  of  the  porosity  profile.  It  was  shown  that 

a porosity  averaged  over  the  porous  layer  depth 
increases  while  Rq  decreases.  This  result  is  rea- 
sonable since  the  quantum  confinement  effect 
tends  to  prevent  dissolution  on  the  local  scale. 

The  fractal  dimension  was  calculated  for 
the  crystal  part  of  the  porous  structures,  formed 
by  computer  simulation  on  three-dimensional 
mesh.  The  dependence  of  fractal  dimension  upon 
size  is  shown  in  Fig.4.  It  changes  from  0.1  to  3 
as  the  size  of  the  structure  varies  from  1.3  to  4 
mesh  steps,  which  corresponds  to  real  sizes  from 
0.3  to  1 nm.  The  X-ray  scattering  experiment 
[18]  showed  that  mesoporous  structures  are 
characterized  by  fractal  dimension  at  sizes  up  to 
100  nm.  Nanoporous  structures  should  display 

the  fractal  properties  at  smaller  sizes. 

3.2.  THE  RESULTS  OF  SIMULATION  OF 
SELICON-ON-POROUS  SILICON  EPITAXY. 

Three-dimensional  images  of  grown  lay- 
ers were  obtained  for  various  porosity,  pore  size 
and  molecular  beam  density,  F.  The  typical 
fragments  of  layer  surface  before,  during  and 
after  Si  deposition  on  PS  are  shown  in  Fig.  5 
(F=l  bilayer/s)  and  Fig.  6 (F=10  bilayer/s).  The 
cross  sections  at  the  edges  with  the  insets  illus- 
trate the  pore  morphology  evolution.  The  differ- 
ence in  the  kinetics  of  epitaxy  for  the  cases  of  various  F is  clearly  seen.  Under  low  flow 
density,  pores  are  the  effective  traps  for  adatoms  and  the  Si  island  formation  is  sup- 
pressed until  pores  are  covered.  At  high  flow  density,  the  adatoms  form  the  islands 
between  the  pores.  The  islands  are  formed  as  far  from  the  pores  as  a diffusion  length. 
This  is  marked  in  case  of  large  distance  between  the  pores.  Thus,  the  porosity  chosen 
for  Fig.  1 and  Fig.  2 is  rather  low.  It  is  an  interesting  result  to  find  that  a change  in 
pore  structure  occurs  only  in  the  upper  (of  1-2  bilayer  thick)  layer. 

In  order  to  specify  the  detailed  mechanism  of  the  Si  growth  over  pores,  the  sta- 
tistical analysis  of  random  atomic  bond  configurations  for  adatoms  inside  pores  was 
carried  out.  Figs.  6 a,b,c  illustrate  the  idealized  scheme  for  formation  of  Si  pendant 
layer.  Circles  show  the  regular  atom  positions  in  (111)  plane.  Fig.  6 d shows  the  frag- 
ment of  diamond-type  structure  elucidating  the  location  of  the  positions  over  depth:  the 


Fig.5.  Fragment  of  PS  (1 1 1)  at 
different  stages  of  Si  epitaxy 
(P= 0.18,  r=1000  K,  F=  1 bilayer/s). 
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large  circles  correspond  to 
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the  upper  surface  layer.  The  positions  occupied  by  atoms  are 
marked  in  black,  those  most  probable  to  be  occupied  - 
in  grey  in  Figs.  6a,b,c.  A white  circle  inside  a black 
one  denotes  an  atom  in  the  inner  bilayer.  At  first,  the 
metastable  nucleus  (a  chain  of  atoms  or  a single  atom 
linked  with  the  structure  by  two  bonds)  is  formed  at 
the  edge  of  a pore  (Fig.  6a).  Under  certain  conditions, 
depending  on  temperature  and  molecular  beam  den- 
sity, the  metastable  nucleus  exists  long  enough  to 
form  a bond  with  a new  atom  arriving  with  the  sur- 
face diffusion  flow.  The  most  probable  positions  to  be 
occupied  (grey  circles  in  the  Figs.  6a,b,c)  are  linked 
with  the  structure  by  two  bonds  at  least.  At  all  stages 
of  pendant  layer  formation  these  positions  appear  to 
be  located  inside  two  upper  bilayers.  Obviously  this 
factor  determines  the  thickness  of  the  pendant  layer 
growing  by  the  perimeter  over  pore.  The  analogous 
series  of  enlarged  fragments  near  a pore  obtained 
during  the  computational  experiment  showed  close 
resemblance  with  those  presented  at  Figs.  6.  The 
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Fig.6.  Formation  of  a pendant  layer.  Occupied  positions  are 
marked  in  black,  those  most  probable  to  be  occupied  next  are 
marked  in  grey.  Figs,  a-c  illustrate  the  step-by-step  occupa- 
tion of  regular  positions  inside  a pore;  d - fragment  of  dia- 
mond type  structure  elucidating  the  location  of  the  positions 
over  depth:  large  circles  correspond  to  the  topmost  surface 
layer. 


similar  double-linked  atomic  configurations  play  the 
role  of  metastable  nuclei  in  the  simulation  process. 

However,  the  “coast  line” 
at  the  edge  of  the  simu- 
lated pores  is  more  jagged 
in  comparison  with  those 
in  Figs.  6. 

The  film  growth 
kinetics  and  changes  in 
overgrown  pore  moiphol- 
ogy  on  atomic  scale  were 
studied.  We  have  intro- 
duced the  effective  depos- 
ited layer  thickness,  w, 
which  is  the  amount  of  Si 
deposited  for  pores  to  be 
overgrown.  The  1/w  value 
characterizes  the  effec- 


tiveness of  overgrowth. 

The  w vs  F is  shown  in  Fig.  7.  The  porosity  is  0.22.  In  the  range  of  high  flow 
densities  (F>2  bilayer/s)  w increases  while  F increases.  The  explanation  is  as  follows. 
Under  high  F the  islands  nucleate  before  the  pores  are  overgrown.  The  higher  F leads 
to  a dense  distribution  of  islands  between  pores.  The  presence  of  islands  decreases  the 
effectiveness  of  the  overgrowth  process  since  adatoms,  which  are  captured  by  islands, 
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Fig.  7.  The  effective  thickness  of  layer  deposited  vs  flow 
density  (P= 0.22,  T=1000  K). 


cannot  form  the  pendant 
layer.  In  the  range  of  moder- 
ate flow  densities  (0.2  bi- 
lay  er/s<F<l  bilayer/s)  w does 
not  depend  upon  F.  Under 
this  condition,  the  rate  of  all 
processes  is  proportional  to  F. 
Under  low  flow  densities 
(F<0. 1 bilayer/s)  w increases 
while  F decreases.  This  effect 
can  be  interpreted  in  terms  of 
the  detailed  mechanism  of 
pendant  layer  formation  de- 
scribed above.  Under  low  flow 
density  the  probability  of  the 
transition  of  the  nucleus  from 
a metastable  to  a stable  state 
decreases,  and  most  of  ada- 
toms migrate  inside  the  po- 
rous layer.  This  decreases  the  effectiveness  of  overgrowth. 

The  atomic  force  microscopy  data  [19]  give  an  evidence  that  the  surface  be- 
tween pores  is  rough  (with  an  average  relief  height  ~0. 1-0.2  nm).  The  initial  smooth 
3D  grid  (160x160x16,  P=0.14)  was  exposed  to  an  artificial  roughening  procedure, 

during  which  the  surface 
atoms  one  by  one  were 
randomly  removed  from 
the  grid.  The  meshes  ob- 
tained at  various  intervals 
were  used  as  the  initial 
substrates  in  simulation  of 
epitaxy.  For  numerical 
characterization  of  porous 
layer  surface  morphology 
the  roughness  function,  S 
[12],  was  calculated  and 
defined  as  the  relative 
number  of  atoms  situated 
in  the  vertical  sections  of 
surface. 

The  dependence  of 
roughness  upon  the 
amount  of  Si  deposited  is 
shown  in  Fig.8  for  various 
initial  surface  morpholo- 


Amount  of  Si  deposited  (bilayer) 

Fig.  8.  The  dependence  of  surface  roughness  upon  the 
amount  of  Si  deposited  for  various  initial  surface  mor- 
phologies (P=0.14,  1=1000  K,  F=  5 bilayer/s). 
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gies  (curves  2-5).  For  comparison,  curve  1 is  plotted  corresponding  to  the  initial  sur- 
face which  is  smooth  between  pores.  At  the  initial  stage,  curves  2-5  have  a sharply 
descending  section,  illustrating  the  effect  of  polishing.  This  effect  arises  due  to  the 
intensive  capture  of  adatoms  in  micropores.  The  next  stage  of  growth  is  defined  by  the 
competition  between  the  formation  of  a pendant  layer  and  the  growth  of  islands  be- 
tween pores.  Arrows  show  the  points  in  the  curves  which  correspond  to  the  overgrowth 
of  pores.  The  interesting  result  is  that  w increases  non-monotonously  with  the  initial 
roughness  S.  The  brief  explanation  of  this  is  as  follows.  Roughening  changes  the  ratio 
between  the  surface  areas  occupied  by  the  upper  and  the  lower  monolayers  compiling 
the  (1 11)  bilayer.  This  leads  to  a critical  change  in  the  effectiveness  of  the  nucleation 
of  islands  and  accounts  for  the  non-monotonous  dependence  of  w upon  initial  S.  The 
dependence  is  not  periodical  because  the  removal  of  monolayers  is  coupled  with  the 
increase  in  surface  roughness. 


4.  Conclusion 

The  present  model  provides  a more  accurate  description  of  PS  formation  by 
electrochemical  etching  in  HF.  The  simulated  structures  exhibit  close  resemblance  to 
those  obtained  in  experiment.  The  3D  meshes,  corresponding  qualitatively  to  nanopo- 
rous  Si  structures,  are  taken  for  simulation  of  Silicon-on-Porous  Silicon  Epitaxy.  The 
simulation  is  carried  out  on  the  basis  of  the  Gilmer  model  extended  for  the  case  of  a 
relief  surface.  The  possible  mechanism  of  overgrowth  of  the  porous  layer  is  proposed 
via  the  pendant  layer  formation.  In  terms  of  this  mechanism  the  dependence  of  growth 
kinetics  upon  the  porosity,  molecular  flow  density  and  initial  surface  roughness  are 
interpreted. 

This  work  was  supported  in  part  by  the  Russian  National  Programs  “Surface 
atomic  structures”  under  grant  97-3.21  and  “Solid  nanostructure  physics”  under  grant 
99-1131. 
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